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Abstract 

Background: Teleost fish present a high diversity of sex determination systems, with possible frequent evolutionary 
turnover of sex chromosomes and sex-determining genes. In order to identify genes involved in male sex 
determination and differentiation in the platyfish Xiphophorus maculatus, bacterial artificial chromosome contigs 
from the sex-determining region differentiating the Yfrom the X chromosome have been assembled and analyzed. 

Results: A novel three-copy gene called teximY (for testis-expressed in Xiphophorus maculatus on the Y) was 
identified on the Y but not on the X chromosome. A highly related sequence called teximl, probably at the origin 
of the Y-linked genes, as well as three more divergent texim genes were detected in (pseudo)autosomal regions of 
the platyfish genome. Texim genes, for which no functional data are available so far in any organism, encode 
predicted esterases/lipases with a SGNH hydrolase domain. Texim proteins are related to proteins from very 
different origins, including proteins encoded by animal CR1 retrotransposons, animal platelet-activating factor 
acetyl hydrolases (PAFah) and bacterial hydrolases. Texim gene distribution is patchy in animals. Texim sequences 
were detected in several fish species including killifish, medaka, pufferfish, sea bass, cod and gar, but not in 
zebrafish. Texim-Wke genes are also present in Oikopleura (urochordate), Amphioxus (cephalochordate) and sea 
urchin (echinoderm) but absent from mammals and other tetrapods. Interestingly, texim genes are associated with 
a Helitron transposon in different fish species but not in urochordates, cephalochordates and echinoderms, 
suggesting capture and mobilization of an ancestral texim gene in the bony fish lineage. RT-qPCR analyses showed 
that Y-linked teximY genes are preferentially expressed in testis, with expression at late stages of spermatogenesis 
(late spermatids and spermatozeugmata). 

Conclusions: These observations suggest either that TeximY proteins play a role in Helitron transposition in the 
male germ line in fish, or that texim genes are spermatogenesis genes mobilized and spread by transposable 
elements in fish genomes. 
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Background 

In contrast to the situation observed in mammals and 
birds, where sex determination systems and master sex- 
determining genes have been conserved over long pe- 
riods of evolution, sex determination is hypervariable in 
teleost fish [1-4]. Related fish species and even popula- 
tions can have different genetic sex-determining systems, 
with various impact of environmental factors. Fish sex 
chromosomes are generally poorly differentiated. They 
are therefore considered to be evolutionary young and 
might have emerged independently in different fish line- 
ages. Hence, parallel studies on different fish models are 
necessary to better understand the evolutionary dynamics 
of sex determination in fish. 

Such studies have demonstrated that different, even 
closely related fish species can have different master sex- 
determining genes. The first master sex-determining gene 
in fish has been identified in the medaka Oryzias latipes. 
This gene, dmrtlbY, is a Y-chromosomal duplicate of the 
autosomal gene dmrtl, which encodes a transcription fac- 
tor with DM domain involved in male development in 
vertebrates [5,6]. Interestingly, dmrtlbY is not present in 
related species from the same genus, which therefore 
must possess other master sex-determining genes. This is 
the case in Oryzias luzonensis, where the putative master 
sex-determining gene is gsdf (gonadal soma-derived 
growth factor), another gene from the sex-determining 
cascade [7]. 

More distant species have also different master sex- 
determining genes. In the Patagonian pejerrey Odontesthes 
hatcheri, a Y-linked duplicate of the amh gene, which en- 
codes the anti-Mullerian hormone belonging to the TGF- 
|3 superfamily, might drive sex determination [8]. In the 
rainbow trout Oncorhynchus mykiss, the master sex- 
determining sdY is a Y-specific duplicate of the interferon 
regulatory factor 9 gene [9]. Added to the results obtained 
in O. latipes, these observations underline the importance 
of gene duplications in the formation of new master sex- 
determining genes in fish. However, the situation is differ- 
ent in pufferfishes from the Takifugu genus. In Takifugu, 
phenotypic sex might be controlled by two alleles of the 
anti-Mullerian hormone receptor type II gene amhr2, 
which is located on both the X and Y chromosomes. A 
single nucleotide polymorphism differentiating both al- 
leles affects one amino-acid in the kinase domain of 
Amhr2. Females are always homozygous for one allele, 
while males are heterozygous [10]. Hence, hyper variability 
of sex determination might be associated with frequent 
evolutionary switch of master sex-determining genes at 
the top of more conserved downstream components of 
the sex determination and differentiation cascade, follow- 
ing the principle of "masters change, slaves remain" [11], 
Further studies on different fish models are necessary to 
test the extent of this hypothesis. 



Another classical fish model to study sex determination 
is the platyfish Xiphophorus maculatus. This species pos- 
sesses three different types of sex chromosomes: X, Y and 
W. Males can be XY or YY, while XX, WX and WY geno- 
types lead to female development [12]. Using a strain 
with XY male heterogametic sex determination, the pos- 
itional cloning of sex-linked genes has been initiated. A 
bacterial artificial chromosome (BAC) genomic library of 
XY males has been constructed, and BAC contigs cover- 
ing the region differentiating the X from the Y chromo- 
some have been assembled and sequenced [13,14], Very 
recently, the genome of an XX female has been se- 
quenced [15]. The sex-determining region was found to 
be prone to DNA rearrangements including deletions, 
transpositions and duplications [12,16]. We report here 
the identification of a novel gene called teximY, which is 
found in there copies on the platyfish Y chromosome but 
absent from the X chromosome. This gene is preferen- 
tially expressed in testis, and has been captured and mo- 
bilized by a Helitron transposon in fish. 

Results 

Three teximY genes are clustered in the sex-determining 
region of the platyfish Y chromosome but are absent 
from the X chromosome 

Four overlapping Y-linked BAC clones (B14, B29, B17 and 
N20; Figure 1) from the Rio Jamapa platyfish BAC genome 
library [13] were sequenced to completion and assembled, 
resulting in a total sequence of 585,694 bp in length. BAC 
clone B29 contained a pseudogene called Y-cript PCR 
analyses demonstrated the presence of Y-cript in all males 
tested (> 200 individuals) but its absence from all females 
(> 200 individuals) (data not shown). Hence, the four as- 
sembled BAC clones are closely linked to the master sex- 
determining gene on the Y chromosome, or may even 
contain it. 

Further sequence analysis revealed the presence of 
three copies of a new gene candidate that was called 
teximY (for testis-expressed in Xiphophorus maculatus 
on the Y) in a region of approximately 300 kb (Figures 1 
and 2). These three genes were detected neither in BAC 
clones specific of the X chromosome nor in the genome 
sequence of an XX female ([15]; http://www.ensembl. 
org/Xiphophorus_maculatus/Info/Index). The presence 
of teximYa, teximYb and teximYc was tested by PCR in 
34 platyfish individuals from the Rio Jamapa population 
(17 males and 17 females). All males tested were posi- 
tive for the three genes, while all females were negative. 
This confirmed the presence of teximYa, teximYb and 
teximYc on the Y chromosome of the platyfish in the 
Rio Jamapa population. TeximY genes are located at dis- 
tances between 80 and 280 kb near Y-cript and linked 
to two copies of a type 4 melanocortin receptor gene 
mc4r [14]. 
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Figure 1 Structure of the teximY gene cluster in the sex-determining region of the Y chromosome of Xiphophorus maculatus. Red 

arrowheads show the position and predicted transcription orientation of teximY copies. The four sequenced BAC clones (B14, B29, B17 and N20) 
are shown as black lines. The distance between teximY copies is given. Locations of the Y-specific Y-cript pseudogene and of two copies of the 
type 4 melanocortin receptor gene mc4r [14] are indicated. 



All three teximY copies were in the same transcrip- 
tional orientation (Figure 1). They showed 95-98% nu- 
cleotide identity at the cDNA level and 89-97% amino 
acid similarity, suggesting recent duplication events. 
Software-based gene structure prediction as well as se- 
quence comparison with platyfish expressed sequence 
tags and sequencing of RT-PCR and RACE-PCR prod- 
ucts indicated that the teximYa and teximYb genes are 
both composed of 15 exons, while teximYc is a truncated 
version lacking the last two 3' exons (Figure 2). Sequence 
comparison of teximYa and teximYb flanking regions re- 
vealed extensive sequence identity for more than 5 kb 
upstream and more than 7 kb downstream of the genes, 
suggesting a duplication of at least 25 kb. This analysis 
could not be performed for teximYc neighboring regions 
due to the missing adjacent genomic sequences. 

BLAST analysis of the platyfish female genome sequence 
[15] revealed the presence of several teximY-r elated 
sequences. These sequences were not present on sex- 
chromosomal BAC contigs covering the region differenti- 
ating the X from the Y chromosome [13,14], They are 
therefore located either on autosomes, or alternatively in 
homologous regions of the X and the Y chromosomes 
(pseudoautosomal regions). One copy called teximl pre- 
sented a high level of nucleotide identity (> 95%) with the 
Y-chromosomal teximY genes as well as a very similar 
gene structure (Figure 2). Hence, teximl might corres- 
pond to the molecular progenitor of teximY sequences, 
or teximl is a (pseudo) autosomal duplicate of one of the 
teximY genes. Three more divergent texim genes were 
also detected in the genome, with one of them, texim2, 
mapping to autosomal linkage group 18. Altogether, four 
copies of texim are present in platyfish (pseudo) autosomal 
parts of the genome, with three additional Y-specific copies 
only in males. 



Texim genes encode putative SGNH hydrolases found in 
fish and other metazoans 

In silico analysis of conserved domains identified a 
SGNH hydrolase domain in the predicted Texim pro- 
teins. This domain, which is found in certain esterases 
and lipases with broad substrate specificity and regiospe- 
cificity [17], is characterized by four invariant catalytic 
residues: the catalytic serine (S) at the N-terminus in 
conserved block I, a glycine (G) in conserved block II, an 
asparagine (N) in conserved block III and the catalytic 
histidine (H) in conserved block IV (Figure 3). An aspar- 
tate residue (D) is often present in block IV. The four 
major residues (S, G, N and H) are important for activity, 
with a different catalytic mechanism when compared to 
common alpha/beta-hydrolases [18]. 

Sequence database analysis identified texim sequences 
in different fish species including killifish, medaka (2 very 
similar copies on chromosomes 10 and 21), pufferfish 
(one copy in Sphoeroides nephelus), sea bass, Atlantic cod 
(one copy) and gar (one copy). No sequence was detected 
in zebrafish. Texim sequences were also found in 
Oikopleura (urochordate, one copy), in Amphioxus (ceph- 
alochordate, three copies) and sea urchins (echinoderms, 
two copies). Texim genes were absent from mammals and 
other tetrapods, suggesting loss of these sequences after 
divergence from the fish lineage 400-450 million years 
ago, or alternatively introduction of texim into the fish 
lineage through horizontal gene transfer. 

All identified Texim sequences form a phylogenetic 
group distinct from other proteins with SGNH domains 
(Figure 4). Within the SGNH hydrolase family, Texim 
proteins are related to proteins from very different ori- 
gins, including bacterial hydrolases, animal platelet- 
activating factor acetylhydrolases (PAFah) and proteins 
encoded by animal CR1 retrotransposons. No reverse 
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Figure 2 Exon/intron structure of texim genes in the platyfish genome. The SGNH hydrolase domain(-coding region) is marked in light blue. 
Introns are represented by broken lines. Predicted start codons are in green and stop codons in red. cDNA and predicted protein sizes are 
provided. Helitron transposon positions are indicated. Question marks indicate lack of information due to absence of sequence data. Both coding 
and protein sequences of texim4 are partial. 



transcriptase was found to be associated with any texim 
sequence, suggesting that texim is not part of a CR1 
retrotransposon. 

Texim is associated with Helitron transposons in fish but 
not in other animals 

For all platyfish texim genes for which enough flanking 
sequence data was available for analysis, a Helitron 



transposon coding region was identified at the close 
vicinity 3' from texim genes (3 kb downstream from 
the stop codon of texim2) or even in the last intron for 
teximY/1 (Figures 2 and 5). Helitrons are eukaryotic 
rolling-circle DNA transposons [19,20]. They are able 
to capture and mobilize gene sequences in certain or- 
ganisms and might therefore play a role in gene evolu- 
tion through exon shuffling and gene duplication 
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Figure 3 Conserved amino-acid blocks in the SGNH domain of predicted Texim proteins. Critical residues important for catalytic activity 
and oxyanion hole residues are indicated, as well as minimum-maximum amino-acid distances between the domains. 



[21,22]. Texim was also found to be associated with 
Helitrons in other fish species including medaka, puf- 
ferfish and sea bass (Figure 5). Association between 
texim sequences and Helitrons transposons was found 
neither in Oikopleura (which does not have Helitrons) 
nor in Amphioxus and sea urchins (which contain 
Helitrons), These observations suggested capture and 
mobilization of an ancestral texim gene by a Helitron 
element in the fish lineage after divergence from uro- 
chordates/cephalochordates. "Free" Helitron transpo- 
sons not associated with texim genes are present in 
zebrafish (which is devoid of texim genes) as well as in 
platyfish and other fish species possessing texim-Heli- 
tron associations. 

TeximY genes are preferentially expressed in testis in 
adult platyfish 

Expression of X. maculatus texim genes was analyzed by 
RT-qPCR in the Rio Jamapa population used to construct 
the BAC library (XY males and XX females, Figure 6), as 
well as in the Rio Usumacinta population (YY males and 
WY females) (Figure 7). The (pseudo) autosomal copies of 
texim, namely teximl, texim2, texim3 and texim4, were 
all detected by PCR in genomic DNA of both males and 
females of the Rio Usumacinta population. From the 
three Y-linked copies, only teximYb was found in Rio 



Usumacinta. In this population, teximYb was detected by 
PCR in both males and females, suggesting localization 
either on the Y chromosome-specific segment or on 
(pseudo) autosomal regions. 

In both populations, teximYa (when present) and texi- 
mYb were both preferentially expressed in testis, with 
teximYb being almost exclusively expressed in the male 
gonad (Figures 6 and 7). We could not detect any ex- 
pression for teximYc, suggesting that this truncated 
sequence corresponds to a pseudogene. Teximl was 
strongly expressed in testis and also in gills in both platy- 
fish populations. In addition, a significant expression was 
found in ovary in the Rio Jamapa population. Expression 
of other texim genes was generally strong in the testis, 
with exception of texim2 and texim3 in the Rio Jamapa 
population, which are predominantly expressed in the 
brain (Figure 6). 

Texim is expressed in spermatozeugmata and late 
spermatids in platyfish testis 

Expression of texim gene in platyfish adult testis was fur- 
ther analyzed by in situ RNA hybridization. Histological 
analysis of the platyfish testis revealed the presence of 
many spermatozeugmata, i.e. the aggregation of mature 
sperm (Figure 8A). Hence, the testis analyzed was at stage 
VI of maturation (according to [23]). The vasa gene, 
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Figure 4 Maximum likelihood phylogenetic tree of Texim-related SGNH hydrolase protein domains found in different organisms. The 

Maximum Likelihood phylogenetic tree is based on an alignment of 1 12 amino-acids from the SGNH hydrolase domains. Bootstrap values are 
indicated (1000 repetitions). Accession numbers: Fish Texim: Lipase_GDSL2_Sphoeroides: Sphoeroides nephelus, protein prediction, AF094698.1; 
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encoding an RNA helicase of the DEAD box protein fam- 
ily, was used as a control for gene expression in germ line 
[24]. According to expression patterns observed in other 
fish species, vasa was expressed in spermatogonia and 
primary spermatocytes in the platyfish (Figure 8B/C). In 
situ hybridization for texim 1/teximY (the probe used does 
not allow to differentiate between teximY and teximl) 
showed expression in late spermatids and spermatozeug- 
mata, i.e. at late stages of spermatogenesis (Figure 8D/E). 

Discussion 

Gene duplications play a major role in the formation of 
novel master sex-determining genes in fish [5,7,9] and 
other organisms [25,26]. New master genes appearing at 
the top of the sex determination cascade are generally 
derived from genes already involved in sexual develop- 
ment, as exemplified by the dmrtlbY master gene of the 
medaka [5,6]. However, exceptions may exist, as sug- 
gested for the immune gene paralog sdY recruited as the 
master sex-determining gene in the rainbow trout [9]. 
Frequently, master sex-determining genes in vertebrates 
encode transcription factors (Sry in therians, Dm-W in 



Xenopus laevis, Dmrtl in Gallus gallus, dmrtlbY m Ory- 
zias latipes), but some recent data have demonstrated 
that it is not the general rule, at least in fish {gsdf in 
Oryzias luzonensis, amhY in Odontesthes hatcheri and 
amhr2 in Takifugu rubripes [7,8,10]). Importantly, all 
master sex-determining genes identified so far in fish are 
expressed in the same manner, in somatic cells sur- 
rounding the germ cells of the differentiating gonads. 

In the platyfish Xiphophorus maculatus, an ovovivipar- 
ous teleost species, the master sex-determining gene has 
not been identified so far. Only a few functional genes 
have been described in the sex-determining region, in- 
cluding the multicopy melanocortin receptor gene mc4r, 
which controls the onset of puberty and is present on 
both the X and the Y chromosomes [14,27]. 

In this study, we have identified a new multicopy gene 
called teximY, triplicated on the Y but absent from the X 
chromosome. TeximY genes are located at distances be- 
tween 80 and 280 kb near Y-cript, a Y-specific pseudogene 
linked to the sex-determining locus. The sex-determining 
region of the platyfish is very rich in repetitive DNA [16]. 
Accordingly, the repeat content of the 0.5 Mb region 
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analyzed in this study is as high as 65.2% (unpublished 
results). In addition, the region containing the master 
sex-determining gene is characterized by a high level of 
genomic plasticity, with frequent rearrangements such as 
duplications, deletions and transpositions [12,16]. The 
presence of three copies of texim is consistent with the 
frequent duplications observed in this region. The auto- 
somal texim 1 gene is a possible molecular ancestor of the 
Y-linked cluster. Y-chromosomal texim copies might be 
the result of mobilization by the flanking Helitron trans- 
posons, or might have been embedded in larger segmen- 
tal duplication events. TeximY genes are preferentially 
expressed in the testis. Hence, the situation observed on 
the sex chromosomes of the platyfish is consistent with 
the observation that new sex-biased genes, especially 
male-biased genes, can be formed on sex chromosomes 
through the duplication of autosomal genes [28]. 

There is so far neither biochemical nor biological 
data concerning the functions of TeximY. Identification 
of conserved SGNH hydrolase domains in predicted se- 
quences strongly suggests that Texim proteins work as 
esterases/lipases [17]. In contrast to master sex- 
determining genes described so far in fish, teximY 



genes are not expressed in somatic cells in the adult 
testis. Hence, teximY genes might not correspond to 
platyfish master sex-determining genes. However, add- 
itional expression analyses remain to be performed to 
test if these genes are expressed during embryonic de- 
velopment at the onset of sex determination. 

We have shown that teximY genes are expressed in late 
spermatids and spermatozeugmata, suggesting that these 
genes might function at late stages of spermatogenesis. 
Interestingly, Texim is related to the platelet-activating fac- 
tor (PAF) acetylhydrolase (PAFah) present in different ani- 
mals. In mammals, PAFah is involved in reproduction 
[29-34], but also in other biological processes such as devel- 
opment, inflammation, hemostasis and apoptosis. This en- 
zyme regulates PAF activity via deacetylation [29]. It reduces 
the level of PAF [33], but can also be involved in signaling 
pathways via interaction with other proteins in different tis- 
sues including testis [32,33,35,36]. Recently, a new testis- 
specific protein isoform of Pafahlb2 with possible role in 
spermatogenesis has been identified, the knock-out of which 
causes male infertility in the mouse [32], 

Texim proteins are also related to prokaryotic pro- 
teins. Hence, either prokaryotic or eukaryotic texim-kke 
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Figure 6 RT-qPCR expression analysis of fex/m genes in adult tissues of Xiphophorus maculatus population Rio Jamapa. All qPCR results 
were normalized against the ribosomal protein L7 housekeeping gene rpl7, the value 1 was assigned to the tissue/organ with the highest level of 
expression. M: males (in blue); F: females (in red). TeximYo, teximYb, teximl and texim4 are preferentially expressed in testis, with teximYb being 
almost exclusively expressed in the male gonad. Teximl is also significantly expressed in ovaries. Texim2 and texim3 show predominant expression 
in brain. 



genes share a very ancient common ancestor, or a hori- 
zontal transfer between prokaryotes and eukaryotes was 
involved in the evolution of texim Aike sequences. In 
addition, Texim proteins show homology to proteins with 
esterase domains encoded by animal CR1 retro transpo- 
sons. These enzymes might play a role in interactions 
with cellular membranes, helping the element in pene- 
trating host cells and possibly facilitating horizontal 
transfer [37]. Texim genes are not parts of retrotranspo- 
sons, since they possess multiple introns and are not as- 
sociated with reverse transcriptase genes. However, we 
cannot exclude that texim sequences have evolved from 
ancestral CRl-like elements. 

Interestingly, texim sequences in fish are associated 
with another type of transposons, rolling circle DNA 
transposons called Helitrons. Such elements have been 



already identified on the sex chromosomes of the platy- 
fish [38]. Helitrons are able to capture and mobilize resi- 
dent coding sequences particularly in plants, and might 
therefore be involved in the evolution of gene functions 
through exon shuffling and gene duplication [21,22], 
The absence of association between texim and Helitrons 
in urochordates, cephalochordates and echinoderms 
suggests that the capture event took place specifically in 
the fish lineage. Texim genes might facilitate Helitron 
transposition and transmission in a way similar to the 
role proposed for the Texim-related esterase in CR1 ret- 
rotransposons [37]. However, most Helitron transposons 
are not associated with texim in teleosts, indicating that 
mobilization and spread of Helitrons do not absolutely 
require Texim even in fish. Patchy distribution of texim 
in fish and other vertebrates suggests either recurrent 
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Figure 7 RT-qPCR expression analysis of tex/m genes in adult tissues of Xiphophorus maculatus population Rio Usumacinta. All qPCR 
results were normalized against the ribosomal protein L7 housekeeping gene rpl7, the value 1 was assigned to the tissue/organ with the highest 
level of expression. M: males (in blue); F: females (in red). All texim genes are strongly expressed in testis, with teximYb and texim2 being almost 
exclusively expressed in the male gonad. Texim! is also strongly expressed in gills, and texim4 in brain. 



loss of texim (e.g., in zebrafish and tetrapods) or horizon- 
tal transfer of texim, possibly in association with Helitrons. 

Conclusions 

Texim genes in the platyfish uncover a novel multigene 
family encoding proteins with SGNH hydrolase domain, 
for which no functional data is available so far. The Y- 
specific copies teximYa and teximYb have no counter- 
parts on the X chromosome and are specifically 
expressed in late spermatids and spermatozeugmata. 
Texim genes are associated with Helitrons in fish ge- 
nomes, suggesting that they have been captured and 
mobilized by these transposons. Two non-exclusive hy- 
potheses can explain these results. On the one hand, 



texim genes might have been maintained in Helitrons 
because Texim proteins facilitate the transposition or 
transmission of these elements. Alternately, texim genes 
might correspond to fertility genes with advantageous 
male functions mobilized and spread by Helitrons trans- 
posons. Such genes, which may have deleterious effects 
on females, are frequently located on the male-specific 
Y chromosome in animals. Fish present a very important 
evolutionary turn-over of sex chromosomes. Transposon- 
mediated mobility of fertility genes might allow the rapid 
fixation of spermatogenesis genes on neo-Y chromo- 
somes, as suggested by the presence of three teximY 
genes on the Y but not on the X chromosome of the 
platyfish. 
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Figure 8 In situ RNA hybridization analysis of vasa and texim 1/Y expression in platyfish adult testis. A. Transversal testis section stained 
with Masson's trichrome. B. Vasa sense probe, C. Vasa antisense probe, D. Texim sense probe, E. Texim antisense probe. SG, spermatogonia; SC, 
spermatocytes; ST, spermatids; SZ, spermatozeugmata. The probe used does not allow to differentiate between teximY and texim 1. Scale bars: 
50 Mm. 

k J 



Methods 

Fish 

Xiphophorus maculatus males and females analyzed in this 
study belong to populations Rio Jamapa (Jpl63A) and Rio 
Usumacinta (Up-2) reared at the Plateau de Recherche 
Experimental de Criblage In vivo (PRECI) of the SFR Bio- 
Sciences Gerland - Lyon Sud (US8/UMS3444, Lyon, France) 
and the fish facility of the Biozentrum at the University of 
Wurzburg (Germany). Fishes were kept under standard con- 
ditions and sacrificed following protocols in accordance with 
regulations from the French Ministry of Agriculture and the 
European Union (agreement number A693870602). 

In silico sequence analysis 

Gene structure and protein sequences were predicted using 
FGENESH with the fish training set (http://softberry.com). 



Intron-exon structure was determined by sequence 
comparison with Expressed Sequence Tags (ESTs) down- 
loaded from Genbank as well as by sequencing of RT- 
PCR and RACE-PCR products. BLAST analyses were 
mainly performed against sequence databases accessible 
from the NCBI (http://www.ncbi.nlm.nih.gov/BLAST/) 
and ENSEMBL (www.ensembl.org) servers. Sequence 
alignments were done using MUSCLE [39] via SeaView 
[40,41] and BioEdit [42]. Phylogenetic reconstructions 
were performed by Maximum Likelihood using PhyML 
4.0 with 1000 bootstrap repetitions under the LG model, 
with the following additional parameters: optimized in- 
variable sites, BioNJ as a starting tree and Best of NNI & 
SPR tree searching options [41]. The VISTA online soft- 
ware was used to perform alignments with mVISTA plot 
in order to visualize the level of sequence conservation 
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(http://genome.lblgov/vista/mvista/submitshtml). Con- 
served protein domains were detected via NCBI (http:// 
www.ncbi.nlm.nih.gov/Structure/cdd/wrpsbxgi), PFAM 
(http://pfam.sanger.ac.uk/) [43] and annotated with Gene 
Ontology terms (http://www.geneontology.org/) [44]. Xipho- 
phorus autosomal and Y-linked texim sequences have been 
deposited in GenBank under accession numbers 
[KF433067- KF433073]. 

Gene expression analysis 

Total RNA was extracted from brain, eyes, gills, liver, 
spleen, muscle, skin and gonads of X. maculatus males 
and females from the Rio Jamapa and Rio Usumacinta 
populations using the TRIZOL reagent (Molecular Re- 
search Center, Inc.). Reverse transcription was per- 
formed with lug of total RNA using RevertAid First 
Strand Synthesis kit and oligo-dT primers (Fermentas). 
Expression patterns were determined by quantitative 
RT-PCR using gene-specific primer pairs. Primers used 
were GCTGCTTGCAGTTTGGATT and GTTTGATG 
GAGCCATGACA for teximYa, AGCCACGGGATGCA 
GAT and GTTCCTGTCCCTTAACCACAAC for texi- 
mYb, AGCCACGGGATGCAGAT and GTTCCTGTCC 
CTTAACCACAAC for teximl, CCTCTGCGGAAGAG 
TTGAAG and G TACG C A AC A A AG CG TC A A A for 
texim2, GGAGACCTGCTGGACTCTTC and TTCCAC 
CTTCGGACAACATA for texim3, and GACGCTCTG 
AGGCAGCTAA and C ATC ACG GGAGG C A AC AT for 
texim4. Analysis of relative gene expression was done 
using the DeltaDeltaCt method [45]. All qPCR results 
were normalized against the ribosomal protein L7 house- 
keeping gene rpl7 ([46,47] Additional file 1: Figure SI). 

In situ RNA hybridization 

In order to generate probes for in situ hybridization, pla- 
tyfish vasa (663 bp, [GenBank: KF424536], primers 
GGTTACCGTGGAAAAGACGA and CACCTTTCCTC 
TCCCAATCA) and teximl (542 bp, primers TTGGCT 
GTCTGTTCTCGTTG and GAGAAACGGACAGGATT 
GGA) partial cDNAs were amplified by PCR from testis. 
Amplified fragments were cloned into the pCRII vector 
(Invitrogen Life Technologies). Digoxigenin (DIG) -la- 
beled antisense RNA probes were synthesized by in vitro 
transcription from 7/mdIII-linearized plasmid constructs 
using the SP6 RNA polymerase, while negative control 
sense RNA probes were produced from X/zoI-linearized 
constructs using the T7 RNA polymerase. In situ hybrid- 
izations were performed on testis sections. Testes were 
fixed in 4% paraformaldehyde (PFA) in Phosphate Buff- 
ered Saline (PBS) buffer at 4°C overnight. After several 
washes in ethanol and then in toluene, testes were em- 
bedded in paraffin. Sections of 8um were cut using a 
microtome and then slides were prepared. After deparaf- 
finization of slides with xylene followed by several 



washes in 100% and 70% ethanol and PBS, testes were in- 
cubated for 20 min with proteinase K (20 ug/mL). Prehy- 
bridization was performed for 1 hour at room temperature 
in the hybridization buffer composed of 50% formamide, 
20X Saline Sodium Citrate (SSC) buffer, 0.1% Tween-20, 
10.5 (ig/ml yeast tRNA, 0.05 mg/ml heparin, 5ug/ml sal- 
mon sperm DNA and IX Denhardts solution in diethyl- 
pyrocarbonate (DEPC)-treated water. Both antisense and 
sense probes were diluted in the hybridization buffer for 3 
min at 80°C. Slides with 100 ul of probe (1 (ig/ml) were 
covered with hybridslips and hybridized in hybridization 
buffer at 70°C overnight. The next day, after removal of 
hybridslips in 5X SSC, slides were washed five times in 
hybridization washing solution (5X SSC, 50% formamide, 
0.1% Tween-20 in distilled water) at 70°C for 30 min, then 
two times in 2X SSC at 70°C for 30 min and four times in 
Tris-Buffered Saline (TBS) buffer for 5 min. Slides were 
subsequently incubated for 1 hour in a blocking solution 
(lmg/ml, Roche), followed by incubation in an anti-DIG al- 
kaline phosphatase conjugate (Roche). The next day, slides 
were washed four times in TBS for 5 min and two times in 
alkaline phosphatase buffer (NTMT; 0.1M Tris-HCl, 0.1M 
NaCl, 1M MgCl 2 , 0.1% Tween-20 in distilled water). Coun- 
terstaining was performed with a BM-purple solution 
(Roche) at room temperature until signal detection. Finally, 
slides were washed in PBS and mounted in Aquatex 
mounting media (Merck Milipore, Germany). Microphoto- 
graphs were taken with Leica light microscope with 20X 
magnification. 

Additional file 



Additional file 1: Figure SI. Basal expression profile of the rp 1 7 gene in 
two populations of Xiphophorus maculatus Rio Jamapa and Rio Usumacinta. 
Two independent sets of cDNA have been generated from adult male and 
female organs: brain, eyes, gills, liver, spleen muscles, skin and gonads. 
Technical triplicates of each cDNA have been performed on the same 
qRT-PCR plate. Experiments were done with the Bio-Rad kit using the 
following PCR program: 40 cycles of 94°C and 59°C. Ct (Threshold cycle) 
values were averaged and standard deviation has been calculated. 
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